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1 The metals in X-ray galaxy clusters 

The X-ray emitting hot plasma in galaxy clusters represents the 80 per cent 
of the total amount of cluster baryons. It is composed by a primordial com- 
ponent polluted from the outputs of the star formation activity taking place 
in galaxies, that are the cold phase accounting for about 10 per cent of the 
cluster barionic budget (see Fig.^). 





Fig. 1. (Left) Updated version of the cluster baryonic pie presented in |1 1 j and 
obtained by using /bar.wMAP = 0.176±0.015 from the best-fit results of the WMAP 
3-yrs data a depletion parameter Y = 0.920 ± 0.023 and the gas and stellar 

mass fractions shown on the right which refer to a sample of 58 clusters with T gas > 4 
keV, redshift in the range 0.14— 1.26 and discussed in Ettori et al. (2006, in prep.). 
The dashed lines indicates the median value of / star = 0.016 and /star //gas = 0.129. 



As reference, a number density of 9.77 x 10~ 2 ions of helium and 4.68 x 
10 -5 ions of iron is expected for each atom of hydrogen in a hot plasma 
with solar abundance (as in £Q; for the most recent estimates in ^U] and 
H3, the number density relative to H is 8.51 x 10~ 2 for He, 3.16 x 10~ 5 and 
2.82 x 10 for Fe, respectively; see Table 1). These values imply a mass 
fraction of (0.707,0.738) for H, (0.274,0.250) for He and (0.019,0.012) for 
heavier elements, accordingly to (pQ, [3]), respectively. 
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I discuss here some speculations, and relevant implications, of the sedi- 
mentation of helium in cluster cores (Sect. |2 details are presented in |14| 'l 
and a history of the metal accumulation in the ICM, with new calculations 
(with respect to the original work in |13|) following the recent evidence of a 
bi-modal distribution of the delay time in SNe la (see Sect. [3J . 



2 Helium in X-ray galaxy clusters 

Diffusion of helium and other metals can occur in the central regions of the 
intracluster plasma under the attractive action of the gravitational poten- 
tial, enhancing their abundances on time scales comparable to the cluster 
age. For a Boltzmann distribution of particles labeled 1 with density n± and 
thermal velocity vth = (2kT ' j 'Axnip) 1 ' 2 in a plasma with temperature kT in 
hydrostatic equilibrium with a NFW potential g(r), the drift velocity of the 
heavier ions 2 with respect to 1 is given by (\22\) 



v se d(r) 
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(1) 



where log A is the Coulomb logarithm. In ^1] , we have studied the effects of 
the sedimentation of helium nuclei on the X-ray properties of galaxy clusters. 
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Fig. 2. (Left) Best joint-fit of the gas density and temperature profile of the de- 
projected data of Centaurus (from 19J with a modified NFW gas profile. (Middle) 
Metal abundance profile with the characteristic central drop observed in spatially 
well resolved cool-core clusters. (Right) H-He (thick dotted line) and H-Fe (thin 
dotted line) equipartition time in the Centaurus cluster. The dots mark the radius 
Tin at which He is accumulated from the regions beyond with a sedimentation time 
represented by the solid line: i sc d = J^ "* dr/v BC d(r). The shaded region ranges 
between a t coo i = 2 x fO 7 yrs and the age of the Universe (details in |14|h 



We have estimated the gravitational acceleration by fitting to the observed 
deprojected gas density and temperature profiles some functional forms that 
well reproduce the central steepening and are obtained under the assump- 
tion that the plasma is in the hydrostatic equilibrium with a NFW potential 
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(left panel in Fig. [21 . The observed gas density and temperature values do 
not allow metals to settle down in cluster cores over timescales shorter than 
few 10 9 years. On the other hand, by assuming that in the same potential 
the gas that is now describing a cool core was initially isothermal, the sedi- 
mentation times are reduced by 1-2 order of magnitude within 0.2r2oo- The 
sedimentation of helium can then take place in cluster cores. Even modest 
enhancement in the helium abundance affects (i) the relative number of elec- 
tron and ions 2* = c M = Y J M l Z l N l = 1 + 2AWV H e + E^H 1 He M i Z i N ^ ( a ) 
the atomic mean molecular weight fi = x ^ 1 ^ z ^ j anc l a n the X-ray 
quantities that depend on these values, such as the emissivity e, the gas mass 
density p g oc /j(1 + cm), the total gravitating mass, M to t oc the gas mass 
fraction, / gas oc p 2 (l + cm)- Moreover, we show in ^1] that if we model a 
super-solar abundance of helium with the solar value, we underestimate the 
metal (iron) abundance and overestimate the model normalization or emis- 
sion measure. For example, with Mjje = 3, the measured iron abundance 
and emission measure are ~ 0.65 and 1.5 times the input values, respec- 
tively. Thus, an underestimated excess of the He abundance might explain 
the drops in metallicity observed in the inner regions of cool-core clusters 
(see, e.g., middle panel in Fig- El - It is worth noticing that the diffusion of 
helium can be suppressed, however, by the action of confinement due to rea- 
sonable magnetic fields (see 0), or limited to the very central region (r < 20 
kpc) by the turbulent motion of the plasma. 



3 Iron in X-ray galaxy clusters 

The iron abundance is nowadays routinely determined in nearby systems 
thanks mostly to the prominence of the K-shell iron line emission at rest- 
frame energies of 6.6-7.0 keV (Fe XXV and Fe XXVI). Furthermore, the few 
X-ray galaxy clusters known at z > 1 have shown well detected Fe line, given 
sufficiently long (> 200 ksec) Chandra and XMM-Newton exposures f |18|. 
[To) . )• It is more difficult to assess the abundance of other prominent 
metals that should appear in an X-ray spectrum at energies (observer rest 
frame) between ~ 0.5 and 10 keV, such as oxygen (O VIII) at (cluster rest 
frame) 0.65 keV, silicon (Si XIV) at 2.0 keV, sulfur (S XVI) at 2.6 keV, and 
nickel (Ni) at 7.8 keV. In ^2|> we infer from observed and modeled SN rates 
the total and relative amount of metals that should be present in the ICM 
both locally and at high redshifts. Through our phenomenological approach, 
we adopt the models of SN rates as a function of redshift that reproduce 
well both the very recent observational determinations of SN rates at z > 0.3 
(ID, 0) and the measurements of the star formation rate derived from UV- 
luminosity densities and IR data sets. We then compare the products of the 
enrichment process to the constraints obtained through X-ray observations 
of galaxy clusters up to z ~ 1.2. 
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Table 1. Adopted values for the average atomic weight (W), solar abundance by 
number with respect to H (A, from pQ and 0; R = ^4aos/^4ag89), total synthesized 
isotopic mass per SN event (mj„ from deflagration model W7 and race integrated 
over the mass range 10 — 50Mq with a Salpeter IMF; see |17p and corresponding 
abundance ratios by numbers with respect to Fe, Yi — rrii / (WiAi) x (WF a AF e )/mF e . 



metal W t!ag89 ^4ao5 R m Ia Y Ia race Ycc 



He 4.002 9.77e-2 8.51e-2 

Fe 55.845 4.68e-5 2.82e-5 

O 15.999 8.51e-4 4.57e-4 

Si 28.086 3.55e-5 3.24e-5 

S 32.065 1.62e-5 1.38e-5 

Ni 58.693 1.78e-6 1.70e-6 



0.87 - 

0.60 0.743 - 0.091 - 

0.54 0.143 0.037 1.805 3.818 

0.91 0.153 0.538 0.122 3.526 

0.85 0.086 0.585 0.041 2.284 

0.95 0.141 4.758 0.006 1.647 



The observed iron mass is obtained from ||] and |2S] (see also ^21) as 

r R 

M F0jObs = 4ttA Fc W Fc / Z Fc (r)p R (r)r 2 dr, (2) 
Jo 

where Z Fc {r) is the radial iron abundance relative to the solar value A Fo , W Fc 
is the iron atomic weight and pn(r) is the hydrogen mass density. To recover 
the history of the metals accumulation in the ICM, we use the models of 
the cosmological rates (in unit of SN number per comoving volume and rest- 
frame year) of Type la, n a , and core-collapse supernovae, rcc, as presented 
in [5] and [231 - We estimate then the iron mass through the equation 

M FeySN = M FcM +M FCyCC = ^2 (™Fc,iana(*) + m Fcfic r cc (dt)) dtdV (3) 

dt.dV 

where mF C ,ia and mFe.cc are quoted in Table 1, dt is the cosmic time elapsed 
in a given redshift range, dV is the cluster volume defined as the volume cor- 
responding to the spherical region that encompasses the cosmic background 
density, p h = 3H%/(8irG) x J? m (k4x 10 10 Af Q Mpc~ 3 for the assumed cos- 
mology), with a cumulative mass of M vir w 6.8(fcT/5keV) 3 / 2 x 10 14 M Q 
We obtain that these SN rates provide on average a total amount of iron that 
is marginally consistent with the value measured in galaxy clusters in the red- 
shift range — 1, and a relative evolution with redshift that is in agreement 
with the observational constraints up to z w 1.2. We predict metals-to-iron 
ratios well in agreement with the X-ray measurements obtained in nearby 
clusters implying that (1) about half of the iron mass and > 75 per cent of 
the nickel mass observed locally are produced by SN la ejecta, (2) the SN 
la contribution to the metal budget decreases steeply with redshift and by 
z w 1 is already less than half of the local amount and (3) a transition in the 
abundance ratios relative to the iron is present at z between ~ 0.5 and 1.4, 
with SN CC products becoming dominant at higher redshifts. 
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Fig. 3. (From left to right, from up to bottom): (A) SN la rates and number 
ratio of SNe CC to la as a function of z for the delay time distribution functions 
4>{td) considered in |S] (dashed lines) and |16l (the solid line is the sum of the 
contribution from "tardy" /squares and "prompt" /dots population). (B) Predictions 
(big dots; using <f>(td) from |16|) of the local metal ratios compared with the X-ray 
measurements in |24| and (C) Accumulation history of the iron mass, Mp e , as 
obtained from eq. (3) with rates from |16|. The solid line shows the sum of SN la 
(dashed line) and SN CC (dotted line) contribution. The accumulation history is 
here compared to the total iron mass measured in samples of local (|3]) and high— 2 
f|25|) clusters for a typical object at 5 keV. (D) Observed (|1]) and predicted (as 
in panel A and normalized to the observed value at z w 0.4) iron abundance as 
function of redshift. (E-F) The last two plots show the metal abundance ratios as 
function of redshift as predicted from a <t>(td) described by a "narrow" Gaussian 
with r = 4 Gyr and a td = 0.2r as in |S] (E; see the corresponding rates in panel 
A) and from a bi-modal SN la distribution as in |16| (F). Horizontal lines indicate 
the redshift region where the SN rates are actually measured. These calculations 
refer to the solar abundance values in £Q for a direct comparison with the data. By 
adopting the compilation in SJ, the trends are confirmed with the most relevant 
change being the predicted Ni/Fe« 2 at z — 0. 

Recently, Mannucci et al. (2006; see also [3U]) argue that the data on 
(i) the evolution of SN la rate with redshift, (ii) the enhancement of the 
SN la rate in radio-loud early type galaxies and (iii) the dependence of the 
SN la rate on the colours of the parent galaxies suggest the existence of 
two populations of progenitors for SN la, one half (dubbed "prompt" ) that 
is expected to occur soon after the stellar birth on time scale of 10 8 years 
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and is able to pollute significantly their environment with iron even at high 
redshift, the other half ("tardy") that has a delay time function consistent 
with an exponential function with a decay time of ~ 3 Gyrs. Due to the 
larger production of SN la at higher redshift (see panel A in Fig. [3J, the 
rates suggested by Mannucci et al. predict a more flat distribution of the Fe 
abundance as function of z than the rates tabulated in |Hj, but still show a 
negative evolution partially consistent with the measurements up to z ss 1.2 
(panel D in Fig.OJ). Furthermore, these rates provide a better agreement than 
the ones in |Kj both with local abundance ratios (panel B in Fig.[3J) and with 
the overall production of Fe (panel C) , that at z — has been released from 
SNe la by more than 60 per cent. A potential discriminant is the behaviour 
of the abundance ratios relative to Fe in the ICM as a function of redshift. 
providing the rates in ^1 values of O/Fe and Ni/Fe that are a factor of 2 
different from the predictions of the single-population delay time distribution 
function (see panels E and F). 
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